Tunneling and Non-Universality in Continuum Percolation Systems 
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The values obtained experimentally for the conductivity critical exponent in numerous percolation 
systems, in which the interparticle conduction is by tunnelling, were found to be in the range of 
to and about to + 10, where to is the universal conductivity exponent. These latter values are 
however considerably smaller than those predicted by the available "one dimensional" -like theory of 
tunneling-percolation. In this letter we show that this long-standing discrepancy can be resolved by 
considering the more realistic "three dimensional" model and the limited proximity to the percolation 
threshold in all the many available experimental studies 

PACS numbers: 72.60.+g, 64.60.Ak 
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While the existence of non-universality in physical 
properties of percolation systems has been established 
some eighteen years ago |l|, |2( the discrepancy between 
the numerous experimental results and the correspond- 
ing available theories is still an unresolved issue |E 0. 
In particular the values of the critical exponent of the 
electrical conductivity, t, reported in hundreds of works 
IE IE IE IE on various composite materials in the 
last twenty years, are not quantitatively accounted for 
by those theories. The major difficulty in the compar- 
ison of the theoretical predictions with the experimen- 
tal results appears to be the lack of the experimental 
geometrical-structural information that is expected to 
yield the diverging (but normalizable) distributions of the 
local-conductances. The latter is the well kno wn El , IE IE 
pre-requisite for the non-universal behavior 10]. The 
essence of the non-universal behavior is that the global 
resultant resistance of a percolation system that is given 
by 



R t oc (p-pcY 



is determined by 



R t oc (R) (p - p c 



-to 



(1) 



(2) 



where (R) is the average value of the "bond" (inter- 
nearest-neighbour conducting particles) resistance in the 
system, p is the bond occupation probability, p c is the 
threshold of the system's electrical connectivity, t is the 
conductivity exponent and to is the "universal" critical 
exponent that is determined solely by this connectivity. 
Since in a random system there is no correlation between 
the geometrical position of a bond and its resistance, any 
random subsystem of p c bonds for a given p, will provide 
a connected-conducting network. In particular, if the 
latter network is chosen by the descending values of the 
bond conductance, g, the value of (R) will be an average 
that is determined by g c , the smallest value of g in that 
sub- network I n the case where the distribution of 

the g values, h(g), diverges as g — » 0, the value of g c 



will diminish as p approaches p c , yielding a diverging be- 
havior of (R). This behavior has been demonstrated by 
Kogut and Straley [l(it (KS) for the distribution: 



(l-a)g- 



(3) 



yielding that (R) = [(1 — a)/a](g~ a — 1). Hence, for the 
non-diverging case (a < 0), (R) is finite while, for (0 < 
a < 1) the diverging (but normalizable) case, one finds 
that |l0|: (R) oc (p — p c )~ tn where, for the distribution 
given in Eq. (JSJ), t n = a/(l—a). Hence the non- universal 
contribution to the conductivity exponent (t n = t — to) 
is determined only by (R) and thus we are able to limit 
our discussion here to the evaluation of this quantity. 

Turning to the problem at hand, i.e. the large quan- 
titative discrepancy between the theoretical predictions 
and the experimental observations, we note that in both 
systems for which theories were advanced, the porous 
media [jj and the tunneling percolation problem [2| (see 
belo w) . the conductors distribution was mapped onto the 
KS [hJ distribution (Eq. ©) yield ing specific predic- 
tions as to the values of t n . However, the experimen- 
tal results were found to be, in general, larger than ex- 
pected in the first class of .systems |2| and smaller or much 
smaller 0, S IE IE IE IE lEl than expected in the second 
class of systems |2l 14 111). One of us has explained 
the first observation by applying phcnomenologically a 
modified h(g) distribution that still retains the KS de- 
pendence (Eq. 0). That approach yielded new limits 
to the theoretical predictions that seem to accommodate 
all relevant experimental data [3 ■ 

In this letter we seek a general understanding of the 
non-universal behavior of the second class of systems, 
i.e. in systems in which the conducting particles are em- 
bedded in an insulating matrix and the transport be- 
tween the particles is by tunnelin g M , \ 4 , y . Ill) . While 
some qualitative explanations [3. l9l llll 1 12) were pro- 
posed for the above mentioned discrepancy in those sys- 
tems, they were unable to account systematically for it 
[3, |ljj • Moreover, only qualitative specific explanations 
have been given [l^, LLZ| to the fact that in some con- 



2 



ducting composites the critical behavior of the transport 
has been found Q, El E| to be composed of a few 
"non-universal" regimes. We note in passing here that 
the dominant network that contributes to R t is that of 
the nearest neighbors' network that yields in practice a 
bona fide percolation system 0] . This is of course 
because of the exponential decay of the tunneling proba- 
bility with the interparticle distance r. Correspondingly, 
the possible divergence of the average "bond" resistance 
(R) is determined by the largest r's of the nearest neigh- 
bours. The corresponding distribution of the nearest 
neighbours distance in the continuum, P(r), that was 
considered previously EI in order to evaluate the ori- 
gin of the nonuniversality of the tunnelling-percolation 
model within the framework of the KS model,was the 
"one dimensional" (ID) Hertz distribution l20j: 



■ exp 



- b 



(4) 



Here r is the distance between the centers of two nearest 
neighbor spherical conducting particle, b is their diam- 
eter and a is the average nearest-neighbour interparti- 
cle (three dimensional) distance that can be estimated 
from 000: ^(a/2) 3 A = 1 where N is the density of 
the conducting particles. Combining this ID distribution 
with the simple (normalized, < g < 1) tunneling de- 
pendence of the interparticle conductance on r 0, E] , 
g(r) = cxp[— (r — b)/d], where d is the typical tunneling 
decay parameter, one obtains an h(q) distribution of the 
form given in Eq. © but with 0,0 tHj: 



a = 1 - 



(5) 



While this ID-like theory, leading to the simple an- 
alytic prediction of Eq.JSJ), appeared to capture the 
essence of the problem and yielded the very simple and 
convenient physical parameter d/(a — b) for the descrip- 
tion of the physical system, it yields t values that were 
much larger than the very many values of t found ex- 
perimentally. For example, the (a — b)/d values in many 
composites 0, 0, 0, 01 are expected El , accord- 
ing to the ID-like phenomenological model, to be of the 
order of 50 (d « 1 nm, a — b w 50 nm) but the highest 
observed 00 values of t in corresponding composites 
was about 10. Another related phenomenon, noticed by 
us, following the compilation of many experimental data 
, is the general trend of the decrease of t with the 
increase of the critical volume fraction of the conducting 
phase, v c , as well as the scatter of the observed t values 
within the to and (about) to + 10 interval for a priori 
similar systems El • O n the theoretical end, while 
the above mentioned simplified model was, for the sake 
of simplicity , a ID- like model, no trials were made to 
test the consequences of this ID-like simplification, both 
conceptually and mathematically, and the findings of one 



of us |2lJ concerning the dependence of h(g) on the di- 
mensionality has not been translated to a prediction of 
t. 

In the present letter we show that the consideration 
of the effect of dimensionality beyond the ID simpli- 
fication yields even a qualitatively different behavior. 
In particular, the corresponding new predictions of the 
higher dimension model enables to account for the above- 
mentioned collection of experimentally observed phenom- 
ena. In fact, we have realized already upon the in- 
troduction of the ID-like model that the probability of 
finding a nearest neighbor should decrease slower than 
the decrease of g with increasing r in order to enable a 
diverging distribution of h(g). This is since if this is not 
the case (i.e. a < d, or a < 0) a non-diverging h(g) 
and thus a universal behaviour of Rt will be obtained. 
In the 3D case the leading term in the P[r) distribution 
[H HI is of the form of exp[-(r 3 - b 3 )/(a - b) 3 } (see 
below) and correspondingly, for any values of (a — b) and 
d, for large enough r, P{r) decreases faster than g(r) 
yielding a non-diverging h(g) as g — > 0, and thus, the 
asymptotic p — > p c critical behaviour is expected to be 
universal. On the other hand this situation cannot be 
mapped onto a simple KS-like result for t n and another 
framework is needed in order to evaluate the t values 
that are to be compared with the experimental observa- 
tions. However, since a encloses the physical information 
of the tunneling-percolation system we keep using it for 
the system characterization in all dimensions. 

The approach we have chosen is the effective medium 
approximation (EM A) ^2^. This choice is justified fol- 
lowing its validity in general 2.3J and for the determina- 
tion of (R) in composite systems |24( in particular, as well 
as our above realization that the contribution t n comes 
only from the average (R). Hence, the fact that the uni- 
versal conductivity EMA exponent, to, is 1, and that p c is 
1/3 in the cubic lattice, rather than the values expected 
from percolation theory (~ 2 and ~ 0.25 respectively), 
simply means that the to value acts as a reference for the 
t n value that we are examining (see below). This yields 
that, considering the tn « 2 value for the universal 3D 
percolation system 0000, the "correct" value of t 
will be larger than the t value derived from our EMA 
calculations just by a unity. 

Turning to the EMA calculation we consider a bond 
percolation model for a cubic lattice with a bond conduc- 
tance distribution function of the form: p(g) = ph(g) + 
(1— p)S(g) where the non-zero conductance values are dis- 
tributed according to h(g). Assuming that the nearest- 
neighbour inter-particle distances r are distributed ac- 
cording to a given distribution function P(r), h(q) re- 
duces to: h(g) = J °°dr P(r)5[g — g(r)], yielding [lfj the 
EMA average bond conductance G as the solution of the 
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FIG. 1: Left panel: Our calculated EMA conductance G as a 
function of p — p c for "dot" particles (b — 0) with a nearest 
neighbour distribution that is given by Eq.Q. The different 
cases refer to different values of the parameter a — 1 — d/(a — 
b). Right panel: the same as for the left panel but for hard- 
core spheres with diameter b > 0. The cases shown have the 
same a values as those in the left panel. 



following integral equation: 
P(r) 



dr 



P-Pc 



g(r) + 2G 2Gp 



(6) 



For a 3D homogeneous dispersion of impenetrable 
spheres of diameter b, P{r) is well approximated by [2fj| : 

24u(7ix 2 +72X + 73) f 3 

P(r) = exp — 8vji(x — 1) 

b 

-I2vj 2 (x 2 - 1) - 24vj 3 (x - 1)] 6{x - 1),(7) 

where x = r/b, < v < 1 is a dimensionless parameter 
(coinciding with the volume fraction of the conducting 
inclusions), 9 is the unit step function and 
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7i 



(i-«) 



3 ' 



72 



v 3 + v 



1 



2(l-w) 3 ' 73 2 (l-v) 3 ' 



(8) 



The parameter v controls the value of the mean nearest- 
neighbor interparticle distance a (that we used in solving 
Eq. ©) through the relation: a = J drrP(r). From 
Eqs. lO and J7J it is clear that the effective medium 
conductance G is governed by the parameter a/b that 
characterizes P(r) and by the tunneling parameter d. A 
numerical solution of the integral equation given in Eq. 
© is plotted in Fig. Q for the case of "dot" particles 
(6 = 0, left panel) and for hard-core spheres with the 
typical d/b = 0.15 value (right panel). The different 
plots of G are shown for different values of the character- 
istic parameter a (see Eq.(j£|). In order to better appre- 
ciate the behaviour of G as p — > p c , the "local" transport 
exponent defined as 



tip) 



rfln(G) 
d\n(p- p c 



(9) 



is plotted in Fig. [21 for the data exhibited in Fig. ^ It is 

clear that for small values of a the local exponent is only 



FIG. 2: The "local" transport exponent t(p), Eq.@, as de- 
rived from the calculated data shown in Fig. Q 



weakly dependent on p and it is very close to t = 1, i.e. 
to the universal value of the EMA. For larger a values, 
t(p) acquires a stronger p — p c dependence which would 
correspond to an apparent non-universality when p is not 
too close to p c . However, as p — > p c , the local exponent 
asymptotically reduces to the universal value of to = 1. 

We have seen that our predicted t values depend on 
the proximity, p — p c , to the percolation threshold p c and 
thus, as we turn to the discussion of the experimental ob- 
servations, we have to compare the parameters that are 
commonly used to quantify this proximity in the contin- 
uum with the above, lattice, p — p c parameter. We note 
in particular that in the latter case p c is of the order of 
unity (0.247 in bond percolation, or 1/3 in EMA, on the 
cubic lattice). On the other hand in the continuum one 
commonly 0, El considers the fractional volume (weight) 
content of the conducting phase, v (w), and its critical 
value v c (w c ). However the latter critical values can be 
vanishingly small 0, H, El EH an d thus, while the abso- 
lute values of v — v c (or w — w c ) may be very small com- 
pared to unity they do not correspond to a close prox- 
imity to v c or w c - Hence, in general, the proximity in 
both cases is better described |2a] in the present context 
by (p — Pc)/Pc and (v — v c )/v c or (w — w c )/w c . Since we 
noted that (p — p c )/Pc ~ (p — Pc), the comparison of the 
t values obtained in the theory and the experiment has 
to be made for the same values of p — p c and (v — v c ) [v c 
or (w — w c )/w c . Moreover, to appreciate the experimen- 
tal "resolution" limits of (v — v c ) or (w — w c ) that are 
achievable thus far in composites, in general, and in sys- 
tems in which the percolation-tunneling model applies 
in particular, let us consider the co-sputtered granular 
metals |4j, |lj, |2jj . For these systems one can typically 
achieve a 50 fold division of the sample with values of v 
that vary from (ideally) to 1. Hence, for the typical 
[27] | v c ~ 0.5 the smallest v — v c interval that can be ex- 
amined, away from v c , is 0.02i> c , yielding that the closest 
proximity of p — p c (p c ~ 0.25, see above) is not smaller 
than about 0.01. As far as we know this is about the 
closest proximity achieved thus far in the study of exper- 
imental systems and thus all the available experimental 
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data correspond to the lattice percolation range of, at 
most, 1 > p ~ p c > 0.01. 

Examining our above EMA results (Fig. [2J in light 
of the above considerations we see that over the above 
widest experimentally achievable p — p c range the value 
of t may be taken as a constant with a deviation of not 
more than ±20 %. Hence, it is not surprising that the 
common fit done in the literature for experimental results 
is taken as representing a single t value over the accessible 
(less than two orders of magnitude) p — p c range, while in 
fact a variable t is present over that range. Indeed, as we 
noted above, indications for the variation of the measured 
transport exponents over the above p — p c range can be 
found in the literature [H0,H>0] Also, the fact that 
the "measured" t values are scattered, and vary between 
very similar composites 0, 0, 0] , but within the limits of 
to and t of the order of to + 10, as suggested here, indicates 
that these observations follow a combination of the small 
variation in internal system parameters (a, 6, d) and the 
limited p — p c intervals that are considered. On the other 
hand, if higher experimental resolutions will be achieved 
(in the preparation of series of samples) a more detailed 
verification of our present EMA predictions, concerning 
the p — p c dependence of t, is expected to be realized. 

As we saw in Fig. |2 the peak in the t values shifts 
to smaller p — p c values with the increase of (a — b) / d, 
i.e. higher t values will be observed the larger the value 
of (a — b)/d (or a), for the accessible p — p c range. For 
this range this is qualitatively similar to the behavior 
to be expected from the ID-like model, but it is by far 
more moderate in this range. However, the most im- 
portant finding is that the t values predicted here are of 
the order of those observed experimentally for the a, b 
and d parameters that characterize the studied compos- 
ites. In fact our present findings that yield relatively low 
(compared to the ID-like prediction to + (a — b)/d — 1) 
t values confirm our above conjecture that it is the, di- 
verging ID and the non-diverging 3D, h(g) distributions 
that are responsible for the different behavior of the ID- 
like and the higher-dimension percolation-tunneling sys- 
tems. The other general trend of the many experimental 
results is that the smaller the v c , (or w c ) values (dimin- 
ishing j2ij even below 0.01) the larger the t value. This is 
well explained now by the above mentioned fact that the 
v — v c values in these composites are much smaller than 
the (proper) lattice-p — p c values, and thus their farther 
proximity to the threshold in these composites, yields, as 
seen in Fig. larger t values. 

In conclusion, we have shown that the values expected 
to be measured for the conductivity exponent, t, in 
tunncling-percolation systems in the continuum, are be- 
tween 2 and the order of 10 for typical ratios of the 
tunneling decay constant and the size of the conducting 
particles. The dependence of t on the proximity to the 
percolation threshold accounts for the many reported ex- 
perimental values and for their scatter between the above 



values, as well as for the increase of the t values with the 
diminishing percolation threshold when it is character- 
ized, as is usually the case, by the fractional volume or 
fractional weight of the conducting phase in the system. 
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